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Breakdown of atomic hyperfine coupling in a deep optical-dipole trap 
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We experimentally study the breakdown of hyperfine coupling for an atom in a deep optical-dipole 
trap. One-color laser spectroscopy is performed at the resonance lines of a single ®^Rb atom for 
a trap wavelength of 1064 nm. Evidence of hyperfine breakdown comes from three observations, 
namely, a nonlinear dependence of the transition frequencies on the trap intensity, a splitting of 
lines which are degenerate for small intensities, and the ability to drive transitions which would be 
forbidden by selection rules in the absence of hyperfine breakdown. From the data, we infer the 
hyperfine interval of the 5Pi/2 state and the scalar and tensor polarizabilities for the 5 P 3/2 state. 

PACS numbers: 42.50.Hz, 32.60.+i, 32.10.Dk, 37.30.-|-i 


Optical dipole traps (ODTs), including optical lat¬ 
tices, are established tools for trapping cold atoms. They 
rely on a position-dependent light shift of the atomic 
ground state. Excited states usually change differently 
so that transition lines are shifted. The line shifts can 
be small, like in shallow traps, or even vanish, like in 
magic-wavelength traps [1, 2]. However, magic wave¬ 
lengths cannot always be employed because either the 
spontaneous emission rate would be too high, or high- 
power lasers are not available, or more than two atomic 
levels are used, making it impossible to find a wavelength 
that is simultaneously magic for all transitions. Never¬ 
theless, the increasing demand for improved control of 
atoms makes it desirable to work in deep ODTs. The 
advantages of deep ODTs include precise localization of 
the trapped atoms, the possibility to perform resolved- 
sideband Raman cooling, and reduced atom loss in the 
presence of heating processes. 

Driving resonant transitions in deep ODTs requires 
precise knowledge of the atomic light shifts. For op¬ 
tical transitions, the excited-state light shifts are non¬ 
trivial and, interestingly, investigations on this subject 
are only at a beginning, although deep ODTs have long 
been employed in the fields of atomic clocks, quantum 
information processing (QIP), and quantum many-body 
physics. Examples include single-atom optical tweez¬ 
ers [3-7], atoms inside high-finesse cavities [8], atoms 
trapped close to nanophotonic waveguides [9] or res¬ 
onators [10], atoms in hollow optical fibers [11], and quan¬ 
tum gas microscopes [12, 13]. Line shifts might also be¬ 
come relevant in future experiments with ions in ODTs 
[14, 15] and molecules in ODTs [16, 17]. 

Here we show that the differential light shifts depend 
nonlinearly on intensity for high intensities, in contrast 
to the well-known linear dependence for low intensities. 
The typical ODT depth, for which the nonlinear effect be¬ 
comes comparable to the natural atomic linewidth is re¬ 
markably small, namely, ks x 0.4 mK for the parameters 
of our experiment, where ks is the Boltzmann constant. 
Moreover, we observe a splitting of resonances which 
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would be degenerate in the linear regime and we observe 
a resonance which would be dipole forbidden in the linear 
regime. The data are obtained by varying the frequency 
of a probe laser near the Di or D 2 line, which illumi¬ 
nates a single ®^Rb atom. When hitting a resonance, 
optical pumping transfers populations between hyperfine 
ground states. After a probe-light pulse, atomic popula¬ 
tions are measured using cavity-enhanced state detection 
[18]. From measured resonance positions at various ODT 
intensities, we infer the hyperfine splitting of the 5Pi/2 
state and the scalar and tensor polarizabilities at 1064 
nm of the 5 P 3/2 state. 

The physical reason for the observed nonlinear behav¬ 
ior is the breakdown of hyperfine coupling between nu¬ 
clear spin and angular momentum of the valence elec¬ 
tron [19, 20]. This is similar to the Paschen-Back effect, 
in which one observes breakdown of hyperfine or even 
fine-structure coupling when applying a strong static 
magnetic field. Hyperfine breakdown has also been ob¬ 
served in a static electric field [21]. Experimental evi¬ 
dence for hyperfine breakdown in a high-intensity light- 
held has been observed in a level-crossing spectroscopy 
experiment [22], which used only one near-resonant light 
held that simultaneously caused and monitored hyperhne 
breakdown. To our knowledge, hyperhne breakdown in 
far-detuned ODTs has not been discussed in the litera¬ 
ture. 

We model the hyperhne breakdown largely analogous 
to the corresponding theory for an electrostatic held [21]. 
We approximate the Hamiltonian as H = Phf + -ffoDT- 
T/hf describes the hyperhne interaction in the absence 
of external helds. It is diagonal in the basis of states 
\J,I,F,mF), abbreviated as [F,m_F). Here, J, I, and 
F are the angular-momentum quantum numbers of the 
valence electron, of the nuclear spin, and of the total 
atom. Their projections onto the z axis are mj, to/, 
and tof- The diagonal matrix elements Ep of Phf in 
the jF, Tof) basis are independent of mp. Throughout 
this work, we use the values measured in Ref. [23] for the 
hyperhne intervals of the 5 F 3/2 state in ®’^Rb. 

Hqbt describes the effect of the ODT light onto the 
atom. The electric held of the ODT light at the atomic 
position is E{t) = Fq cos(a;t) with amplitude Eq and 
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FIG. 1. (Color online) Energy levels of the ?)P 3/2 excited-state 
manifold in ®^Rb in the presence of light shifts created by tt- 
polarized 1064 nm light. Green, blue, purple, and red solid 
lines correspond to \mF\ = 0, 1, 2, and 3. At low intensity, 
E is a good quantum number, whereas at high intensity rtij 
becomes a good quantum number, similar to the hyperhne 
Paschen-Back effect. Dotted lines show linear approximations 
for low and high intensity. In the high-intensity region, the 
set of steeper lines correspond to |»nj| = |. Note that an 
intensity of 1 MW/cm^ corresponds to a ground-state ODT 
depth of ks x 1.5 mK. The \mF\ = 0 and 1 lines of the E = 2 
and E = 3 manifolds are hardly resolved. 


angular frequency lo. We assume that the light is lin¬ 
early polarized, and choose the z axis along Eq (tt po¬ 
larization). Considering the limit where all excited-state 
populations are small, the effect of the light is reduced 
to generating light shifts. In addition, considering the 
limit where all detunings are much larger than all hy- 
perfine splittings, we approximate ffoDT as being diag¬ 
onal in the basis of states |J,/, (abbreviated as 

I'fTT'j,'fTT'i}) with diagonal matrix elements 




^ J(2J-1) 


( 1 ) 


for J > 1 and E^j = —^a^Eg otherwise. This can be 
regarded as analogous to Ref. [21] combined with time 
averaging {E'^{t))t = Eg/2. Note that these matrix el¬ 
ements are independent of mi. Here, and are 
the scalar and tensor polarizabilities and 7 abbreviates 
the set of quantum numbers n, L, and J, where n is 
the principal quantum number and L the orbital angular 
momentum quantum number of the valence electron. An 
additional term containing the vector polarizability [24] 
vanishes because we consider linear polarization. As the 
light is far detuned, the imaginary parts of the polariz¬ 
abilities are negligible. 

In the low-intensity limit, i?HF dominates and i?oDT 
can be treated in first-order perturbation theory. Here, 
H remains diagonal in the jE, mi?) basis. This yields low- 
intensity light shifts of the form — ja.y_F,mF-®o + ^(-®o) 



FIG. 2. (Golor online) Scheme of the experimental setup. 
A single ®^Rb atom is held in a two-dimensional optical lat¬ 
tice inside a Fabry-Perot resonator (gray). The relevant light 
shifts are created by the 1064 nm light. Probe light is applied 
for in-trap spectroscopy. The cavity and an additional laser 
(not shown) are used for hyperhne state detection. 


with the low-intensity polarizability 

s T 3m|. - F{F + 1) 
Q^-y,F,m_F — F p(2F -- 1) 


( 2 ) 


for E > 1 and = o:^ otherwise. Here, 

_ T 3X(X — 1)—4F(F+1)J(J+1) -i-u y — 

®7,F ~ *^7 (2F+3)(2F+2)J(2J-1) With A — F (r + 1) + 

J(J -b 1) — /(/ + 1) [25]. For later reference, we note 
that this implies Oi^pg /2 f =2 = ^ for ®^Rb {I = |). The 
emergence of ot^ p from is somewhat reminiscent of 
the emergence of the hyperhne Lande factor qf from the 
hne structure Lande factor gj. 

For intermediate intensities, the hyperhne coupling 
breaks down and F[ is diagonal in neither the jE, mF) 
nor the \mj,mi) basis. Here, FI can typically be di¬ 
agonalized only numerically and the light shifts depend 
nonlinearly on Eg. For high intensities (but not so high 
that hne-structure breakdown would occur), Hqbt dom¬ 
inates and Hup can be treated in hrst-order perturbation 
theory. Here, H is diagonal in the \mj, mi) basis and the 
energy eigenvalues are linear in Eg up to an offset. 

Figure 1 shows the energy eigenvalues of H for the 
5 E 3/2 manifold of ®^Rb as a function of the intensity 
Iq = ^c£oEq of a plane traveling light wave, where c 
is the vacuum speed of light and eg the vacuum permit¬ 
tivity. This calculation is based on the theoretical pre¬ 
dictions o:^pg /2 ~ —1114(16) a.u. and C(gpg /2 ~ 551(5) 
a.u. [26] for 1064 nm [1 atomic unit = 1.648 78 x 10“"^^ 
J(V/m)“^]. The figure shows the linear regimes for low 
and high intensities as well as the nonlinear behavior 
in the intermediate regime. Note that c/^pg ^2 f =2 “ ® 
causes identical low-intensity behavior for the complete 
E = 2 manifold, whereas different |mF| values within the 
E = 1 or E = 3 manifolds feature different initial slopes. 

Figure 2 shows a scheme of the experimental setup. A 
772 nm light is coupled into the TEMqo mode of a high- 
finesse Fabry-Perot resonator. This light creates a re¬ 
pulsive, one-dimensional (ID) optical lattice for ground- 
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FIG. 3. (Color online) Measured excitation spectrum of a single atom for transitions from the 5S'i/2, F — 1 ground state 
to 5 P 3/2 excited states. A probe light pulse transfers population into the F = 2 ground state by optical pumping, whenever 
resonant with an atomic transition. The final populations are measured. Five atomic resonance lines are clearly resolved. The 
line shows a fit of the sum of five Lorentzians with independent amplitudes, widths, and center frequencies. The boxes label 
energy eigenstates that we assign to the 5 P 3/2 excited states causing these resonances, where terms with amplitudes below 0.01 
are omitted. Without hyperfine breakdown, this frequency range would show only three lines, because the two rightmost lines 
would be degenerate and the leftmost line would be a forbidden transition. 


state ®^Rb atoms with a potential height of As x 0.73 mK. 
Additionally, a 1064 nm light beam crosses the resonator 
center. This beam has a waist (1/e^ radius of intensity) 
of w = 16 ^m and is tt polarized. It is retroreflected and 
creates an attractive, ID lattice. A single ®^Rb atom is 
loaded into the resulting two-dimensional (2D) lattice. 
The atom is located near an antinode of the 1064 nm 
light and near a node of the 772 nm light. Hence, at the 
position of the atom, light shifts from the 772 nm light 
are small (see the Appendix). The relevant light shifts 
for our measurements come from the 1064 nm light. Hop¬ 
ping of the atom between lattice sites is negligible. 

An additional traveling-wave probe light beam with a 
wavelength near 780 nm or 795 nm crosses the cavity 
center perpendicularly to the cavity axis and subtending 
an angle of 45° (out of the image plane in Fig. 2) with 
the 1064 nm light. This light is linearly polarized, per¬ 
pendicularly to the z axis, thus having equal fractions of 
and (j~ polarization. The probe light is frequency 
stabilized to a commercial frequency comb. All cavity 
resonances are far detuned from the probe light, so that 
the cavity does not affect the spectroscopy signal. 

To perform spectroscopy, we first prepare the atom in 
the F = 1 manifold of the 551/2 ground state with optical 
pumping using additional laser beams, not shown in Fig. 
2 . At this point, the atoms are approximately equally 
distributed among the mp substates of this F = 1 man¬ 
ifold. Second, we expose the atom to a 2 /rs long pulse 
of probe light at an intensity of roughly 0.6 mW/cm^ 
and at a fixed frequency. Third, we measure whether 
the atomic population is still in the F = 1 manifold or 
whether the probe light caused a spontaneous Raman 
transition into the F = 2 manifold. This measurement 
uses cavity-enhanced state detection [18]. 

After loading the atom and after each such spec¬ 
troscopy sequence, we perform polarization-gradient 
cooling using additional light beams to cool the atom 
to ^ 30 fiK. During this cooling process, we take yz- 


plane images of the cooling light scattered from the atom 
with a digital camera. With these images, we monitor 
whether the atom is still in place. A complete cycle for 
spectroscopy and cooling takes 0.5 ms. Such cycles are 
repeated for different frequencies of the probe light, until 
the atom is eventually lost. At that point, a new atom is 
loaded. 

Figure 3 shows a measured spectrum. The 1064 nm 
standing-wave ODT has a traveling-wave power of P = 
2.3 W, corresponding to an ODT depth of kp x 3.4 mK. 
The data clearly show five resonance lines. The excited 
states belong to the F = 0, 1, and 2 manifolds and are 
easily found in Fig. 1, in which these data correspond 
to 2.2 MW/cm^. The widths of the resonances are ^ 
1.5 times the natural linewidth of the excited state. We 
attribute most of this broadening to the thermal position 
distribution of the atom (see the Appendix). 

For the right part of Fig. 3, the hyperfine-coupled 
model Eq. (2) predicts three degenerate lines for the 
F — 2 manifold because f =2 ~ appear¬ 

ance of two resolved lines in this region unambiguously 
demonstrates hyperfine breakdown. 

The leftmost line in Fig. 3 belongs to the F = 0 mani¬ 
fold. The F = 0 excited state cannot decay to the F = 2 
ground state because of dipole-selection rules. Hence, 
this line would be invisible to our spectroscopy method 
according to the hyperfine-coupled model. But hyperfine 
breakdown gives the F = 0 excited state an admixture 
of the F = 2 exited state which can decay to the F = 2 
ground state, making the line visible. Note that a similar 
argument applies to all F = 3 states. 

Figure 4 shows a quantitative analysis of the resonance 
positions. The values of the probe detunings at the line 
centers (dots) were extracted by fitting the sum of inde¬ 
pendent Lorentzians to spectra like in Fig. 3 for different 
values of the traveling-wave power P of the standing- 
wave ODT at 1064 nm. Statistical error bars are ~ 0.1 
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FIG. 4. (Color online) Quantitative analysis of differential 
light shifts. Solid lines show fits to the experimentally ob¬ 
served line centers (dots). Dotted lines show linear approxi¬ 
mations for low power. The data cover the (a) 55'i/2 5Pi/2 

and (b) 5Si/2 5 P 3/2 transitions in ®^Rb. Note that unlike 

Fig. 1, these data show differential light shifts. 


MHz, much smaller than the symbol size. We label the 
detuning as zero, where the light is resonant with the 
transition from the F = 1 ground state to the barycenter 
of the excited-state hyperfine manifold in the absence of 
light. 

Figure 4(a) shows data for the 5S'i/2 <->■ 5Pi/2 tran¬ 
sition. In the framework of Eq. (1), the tensor polariz¬ 
abilities for both involved states vanish. Hence, we ex¬ 
pect that the resonance positions depend linearly on Eq 
with a slope of — |Aai/ 2 , independent of F. Here, we 
abbreviated Aaj = a§pj — the experiment, 

there is some uncertainty in the power calibration, in the 
measured waist, and in the geometrical overlap of the in¬ 
coming and retroreflected 1064 nm beams with the atom. 
We express this by replacing Eq tjEq in the theory, 
introducing a dimensionless parameter ry, which would 
ideally be unity. A fit to the data yields the best-fit val¬ 
ues ryAai /2 = —1675(3) a.u. and {Ep ^2 — Ep=\)l2'K‘h = 
816.6(2) MHz. Only statistical uncertainties are quoted 
throughout this text. A thorough analysis of systematic 
uncertainties is beyond the scope of the present work. 
The value for {Ef =2 — EF=i)/‘2,Trh agrees well with the 
measured value 816.66(3) MHz from Ref. [27]. This seems 
to create serious doubt in the measured value 812.29(3) 
MHz from Ref. [28], which deviates from our value by 
~ 20cr. 

Figure 4(b) shows data for the 55'i/2 5 P 3/2 tran¬ 

sition. A fit of the model based on Eq. (1) to the data 
yields the best-fit values r]AaQ /2 = —1590(3) a.u. and 
a^PQl^lAaQ /2 = —0.312(4). Combination with the best- 


fit value from Fig. 4(a) yields Aai/ 2 /Aa 3/2 = 1.054(2). 

We compare these results with the theory values 
« 5 Si /2 = 687.3(5) a.u. [24], afpi /2 = -1226(18) a.u. 

[26], and the above theory values for all for 1064 

nm. They yield the theory predictions AaQ /2 = 

—0.306(4) and Aai/ 2 /^ck 3/2 = 1.062(14), in good agree¬ 
ment with our measured values. 

As long as ry is unknown, our experimental data yield 
only ratios of polarizabilities. The simplest way to pro¬ 
ceed is to assume 7y = 1, which immediately yields ex¬ 
perimental values for Aai/ 2 , Aq; 3 / 2 , and Qf^p 3 / 2 - 
other way to proceed is to use the above theory values 
for ‘^fpi/ 2 ’ which yields Q;fp 3/2 = —1128(18) 

a.u., a^pQf 2 ~ 568(9) a.u., and ry = 0.876(8). The lat¬ 
ter reveals the size of the systematic uncertainty caused 
when one would set ?y = 1. 

Finally, we note that a measurement of the tensor po¬ 
larizability c^pQi 2 does not require the observation of 
hyperfine breakdown, because the light shifts in the low- 
intensity limit already reveal oi^pQf 2 F’ from which one 
can easily calculate (^pQi 2 - 

To conclude, we experimentally observed hyperfine 
breakdown induced by off-resonant light. This splits res¬ 
onances, lifts selection rules, and creates shifts which are 
nonlinear in light intensity. The line shifts are relevant for 
precision experiments and high-fidelity QIP experiments, 
which drive resonant optical transitions in ODTs. The 
lifting of selection rules is relevant for all schemes that 
rely on the forbidden character of a transition, such as 
optical pumping schemes that populate anruF = 0 state. 
The observed effects should be even more dramatic for 
atoms with smaller hyperfine splittings, like lithium. For 
example, if a recent ®Li quantum gas microscope exper¬ 
iment [13] moved from the Di to the D 2 line, it would 
immediately operate deeply in the Paschen-Back regime. 

A related experiment was simultaneously performed at 
the University of Michigan [29]. 
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Appendix 

The spatial distribution of the atoms has a nonzero 
thermal width. As atoms at different positions experi¬ 
ence different light shifts, the resonance lines in our ex- 
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periment are broadened. As the broadening mechanism 
is not symmetric, this also causes a small shift of the line 
center. In the following, we describe a simple model for 
this broadening mechanism. 

We can safely approximate the ground-state trapping 
potential Vg{x) as harmonic, because the atomic tem¬ 
perature is much smaller than the depth of the trap¬ 
ping potential. This yields Vg(x) = Vg(0) -I- Ug{x) with 
Ug{x) = X]i=iwhere m is the atomic mass 
and the uji are the angular frequencies of the ground-state 
trap. Similarly, we approximate the excited-state poten¬ 
tial as Ve{x) = Ve(0) -I- Ueix) with Ue{x) = 
with certain coefficients bi which are typically negative 
in our experiment. 

For a potential depth of kb x 3.4 mK of 1064 nm 
light and a potential height of ks x 0.73 mK of 772 nm 
light, we obtain trapping frequencies of {ujx, ujy,uJz)f2Tr = 
(11,760,480) kHz. For ksT hiOi for all i one could 
safely use a semiclassical approximation to calculate the 
thermal atomic position distribution. The typical atomic 
temperature in our experiment T ~ 30 /iK corresponds 
to ksT ~ 27r7i x 600 kHz so that we are near the border 
of the ksT 3> fioji regime. For simplicity, we neverthe¬ 
less use a semiclassical approximation for the position 
distribution of the single atom 


and A(a;) = Wprobe — ijJiesi'x) the detuning of the probe 
light from resonance, which depends on the position x 
because of the position-dependent differential light shifts 


A(a;) 


Ao-b 


Ugjx) - Ue{x) 

n 


(A.3) 


where Aq = Afree + bg( 0 ) ~ ^e{0) is the detuning at the 
trap center and Afree is the detuning for a free atom, in 
the absence of ODT light. 

Let r 2 denote the partial rate coefficient for sponta¬ 
neous radiative decay of state |e) into the F = 2 ground- 
state manifold. Then averaging over the atomic density 
distribution yields 



For simplicity, we restrict our considerations to a 
regime where the product Iptp of the intensity Ip and 
duration tp of the probe pulse is so low that we can ap¬ 
proximate pgg as time independent. Hence, 9*^2 is also 
time independent and the quantity that we measure in 
our spectra is N 2 {tp) = tpdtN 2 . As we are not partic¬ 
ularly interested in the absolute heights of the spectral 
lines, it suffices to consider the dimensionless spectrum 


n{x)=noe-^^o‘'^\ (A.l) 

where Hq is the peak density and /? = l/fc^T. The 
normalization condition f d^x n(x) = 1 yields no = 
(TOW^/27rA:BT)^/^, where w = is the har¬ 

monic mean of the trapping angular frequencies. Here, 
we assumed that different spatial directions are thermal- 
ized, which we did not verify experimentally. We mea¬ 
sured the temperature only along the x axis. For the 
density distribution (A.l), the standard deviations of the 
atomic position are (a^, cry, Uz) = (750,11,18) nm. Grav¬ 
ity points along the x axis with a gravitational accelera¬ 
tion oi g = 9.8 m/s^. This causes a gravitational sag of 
Aa; = = 2 nm, which is negligible compared to <7x- 

To see how this atomic position distribution affects the 
spectroscopy signal, we need to calculate how the num¬ 
ber N 2 of atoms in the F = 2 ground-state manifold 
changes in time because of optical pumping. Because 
the 2 ps probe pulse is much longer than the natural 
lifetime l/F = 26 ns of the excited state, we can use a 
rate-equation model to describe the temporal evolution 
of the atomic populations. 

If the probe light near-resonantly drives the transition 
from a ground state \g) with F = 1 and a given mp to 
an excited state |e), then the steady-state population in 
internal state |e) is 

Pee{x) = pgg^^ l + [2A(a;)/r]2’ 

where pgg is the internal-state population in |g), A the 
saturation intensity for the |g) -f-)- |e) transition. Ip the 
probe-light intensity (for which we assumed Ip A), 


5(Ao) = 


2L 


Ip^2Pgg 


dtN2 = no d^x 




l + [2A(x)/F]2’ 

(A.5) 


which is normalized to f SdAg = 7 rr/ 2 . 

Using the harmonic approximation, we obtain 

^ ^ 1 + 4[j + (c/g - u,)/hrp ’ 

where we introduced a dimensionless detuning and the 
dimensionless coordinates 



In general, the dependence of the light shifts Vg/e on 
trap intensity is nonlinear. However, for low enough 
atomic temperatures, the thermal atomic distribution 
samples only a small intensity range which is near the 
intensity at the trap center. Hence, we can safely ap¬ 
proximate the light shifts Vg/f^{x) as first-order Taylor 
polynomials in the intensity around the intensity at the 
trap center, yielding 

2 1 

Ug/e{x) = ^ --ag/e,i[£lo.*(a5) - (A.8) 

i=l 

with Eo,i{x) associated with the \i where (Ai,A 2 ) = 
(1064,772) nm. Here we used that the trap intensity 
is proportional to Eq{x). Our notation expresses the 
first-order Taylor coefficient in terms of Og/^ i which has 
the dimension of a polarizability and turns into (ig/e,i = 
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Q^g/e,i for the special case, where the light shift actually 
is linear in intensity. 

For our experiment, the 1064 nm light provides the 
harmonic confinement along x and y, whereas the 772 nm 
light provides the dominant harmonic confinement along 
z. In principle, the 1064 nm light provides additional 
harmonic confinement along z but that is small and we 
neglect it. Hence 

l3{Ug - Ue) = ar,i{ul + ul) + ar,2ul, (A.9) 
where we abbreviated 


^r,i — 1 


(A.IO) 


In cylindrical coordinates, this yields the final result 
2 


5* = 


dUn 


du. 


2 2 
u„e “<? 


1 + 4(^ + TiU? + r2u2)2 ’ 


(A.1I) 


where we introduced two dimensionless temperatures 


n = 


keT 

TiT 




(A.12) 


The integral (A. 11) can be evaluated numerically. 

Table I shows theory values for the scalar and tensor 
polarizabilities relevant for our experiment. The values 
for 1064 nm are the theory values from our paper. We 
calculated the values for 772 nm ourselves in second-order 
perturbation theory, including transitions 55'i/2 <->■ n'Pj' 
and 5Pj e-)- n'5'i/2, and 5Pj <->■ n'Dji using resonance 
frequencies and oscillator strengths from Refs. [30-32]. 

Table II lists the quantity j = 1 — 
which is related to calculated from Table I. Of 
course, one could calculate Eq. (A. 11) numerically for 
each excited state. Here, however, we consider only one 
example, namely, for the typical values ar^i = 2.7 for 
1064 nm and 5^,2 = 1-2 for 772 nm. Combination with 
ksT/hr = 0.10 yields 


Ti = 0.27, T2 = 0.12. (A.13) 


TABLE 1. Polarizabilities in 

a.u. for ®^Rb. 


®5Sl/2 

g 

®5Pl/2 

g 

^5P3/2 

^5P3/2 

1064 nm 687.3(5) 

-1226(18) 

-1114(16) 

551(5) 

772 nm -11995 

1450 

3460 

-4842 


TABLE 11. Values of af// = 1 
on Table 1. 

S/T , s 
'^5Pj/*^5Sl/2 

for ®^Rb based 

A 

^r,l/2 

®r,3/2 

^^3/2 

1064 nm 

2.78 

2.62 

0.20 

772 nm 

1.12 

1.29 

0.60 



Ao/r 


FIG. 5. (Color online) Expected spectral line shapes. The red 
solid line shows the expected line shape numerically calculated 
from Eq. (A. 11) for the typical values for our experiment from 
Eq. (A.13). This line is broadened and shifted because of the 
nonzero width of the thermal atomic position distribution. 
The black dotted line is a Lorentzian fit to this. The green 
dashed line shows the Lorentzian reference for ri = r 2 = 0, 
corresponding to a classical particle located exactly at the 
trap center. 


In Fig. 5, numerical results for the integral (A.11) for 
these values of ti and T 2 are shown as a red solid line. 
The black dotted line shows a fit of a Lorentzian to this 
result. The best-fit values are a full width at half maxi¬ 
mum (FWHM) of 1.25r and a center position of —0.30r. 
For reference, the green dashed line shows the trivial 
Lorentzian obtained for ti = T2 = 0, which has a FWHM 
of F and is centered at 0. The latter can be regarded 
as the zero-temperature limit of Eq. (A. 11), although in 
this limit, the semiclassical approximation is expect to 
be poor. 

Comparison of the measured line widths of ~ 1.5F with 
the width 1.25r from the fit suggests that the line broad¬ 
ening observed in the experiment is to a good part caused 
by the thermal atomic position distribution. 

We emphasize that the line shifts and line broadenings 
in this model are independent of the depth of the 1064 
nm and height of the 772 nm potential, as long as the 
traps are strong enough that the harmonic approximation 
works well and as long as they are weak enough that the 
semiclassical approximation works well. 

The line shift found in this model suggests that all 
lines measured in our experiment are systematically red 
shifted by ^ 0.3r ^ 1.8 MHz relative to the free-space 
resonance lines. Note that within the approximations 
used so far, this has no effect on the polarizabilities and 
hyperfine splittings that we extract from the data. Those 
values will only be affected if additional corrections are 
taken into account, such as the fact that the actual val¬ 
ues of dry are somewhat different for different atomic 
states. Because of the lack of a tensor polarizability in 
the 5 S'i/2 5Pi/ 2 transition, however, this is not the 
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case for values extracted from those data, including the 
5 Pi/ 2 hyperfine splitting. 

To experimentally test the sensitivity to a change in 
the 772 nm potential height, we measured the resonance 
line involving the F = 1771^1 = 1 excited state at two 
potential heights, namely ks x 0.55 and 1.1 mK. The 
relative shift of the line center between these two experi¬ 
ments was measured to be 79(91) kHz, which is consistent 
with zero. 


where [33] 

w{z) = e~^ eTic{—iz) (A.16) 

is the Faddeeva function, erfc(z) = 1 — erf(z) the com¬ 
plementary error function, and erf(z) = ^ dt the 

error function. Equation (A.15) can be derived using 
partial-fraction decomposition and the integral represen¬ 
tation [33] 


Finally, we note that if only one wavelength was used 
for the optical lattice, the integral (A.11) could be solved 
analytically. To this end, one could rewrite it in spherical 
coordinates as 


2iz r°° 

wiz) = — Im2>0. (A.17) 

TT Jo Z — t 

Note that Eq. (A.17) can also be used to analytically 
solve the Uz integral in (A.11), yielding 


5* = 


dUr 


1 -I- 4((5 -I- TiulY 


(A.14) 


A 




dupC '^puplz{up) 


(A.18) 


For (5, Ti € R with ti ^ 0, this integral has the analytic 
solution 


S = '{^zw{z') F C.C., 
2ti 


2; = 


7 - 2(5 

2ti 


sign(ri), (A. 15) 


with 




li-25-2TiV? . 

-Tt2- 

(A.19) 
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